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Smart polymers are a modern class of polymeric materials that often exhibit unpredictable behav¬ 
ior in mixtures of solvents. One such phenomenon is co-non-solvency. Co-non-solvency occurs when 
two (perfectly) miscible and competing good solvents, for a given polymer, are mixed together. As a 
result, the same polymer collapses into a compact globule within intermediate mixing ratios. More 
interestingly, polymer collapses when the solvent quality remains good and even gets increasingly 
better by the addition of the better cosolvent. This is a puzzling phenomenon that is driven by 
strong local concentration fluctuations. Because of the discrete particle based nature of the inter¬ 
actions, Flory-Huggins type mean field arguments become unsuitable. In this work, we extend the 
analysis of the co-non-solvency effect presented earlier [Nature Communications 5, 4882 (2014)]. 

We explain why co-non-solvency is a generic phenomenon that can be understood by the thermo¬ 
dynamic treatment of the competitive displacement of (co)solvent components. This competition 
can result in a polymer collapse upon improvement of the solvent quality. Specific chemical details 
are not required to understand these complex conformational transitions. Therefore, a broad range 
of polymers are expected to exhibit similar reentrant coil-globule-coil transitions in competing good 
solvents. 


I. INTRODUCTION 

The microscopic understanding of smart polymer con¬ 
formations in a mixture of solvents is scientifically chal¬ 
lenging [lj and, at the same time, possesses great techno¬ 
logical implications that span over a broad range of dis¬ 
ciplines [a-lg. Therefore, establishing the links between 
smart polymer conformations and its specific interactions 
is a key to develop any fundamental understanding of 
their solubility. Examples of the most commonly known 
smart polymers include: poly(N-isopropylacrylamide) 
(PNIPAm), poly(N-isopropylmethacrylamide) (PNIP- 
MAm), poly(N,N-diethylacrylamide) (PDEAm), poly(N- 
vinlycaprolactam) (PVCL), and poly(acryloyl-L-proline 
methyl ester) (PAPOMe). When some of these polymers 
are dissolved in a mixture of solvents, such as aqueous al¬ 
cohol solutions, they show a puzzling coil-globule-coil sce¬ 
nario [6pl0, |U|. This interesting phenomenon is termed 
as co-non-solvency. 

Theoretical understanding of these complex phenom¬ 
ena is mostly restricted to a limited number of com¬ 
puter simulation studies Q], 0, LUHlJk which usually deal 
with chemically specific details [alii, [HI 13 • More¬ 
over, these simulations require careful parameterization 
of force fields that can be cumbersome, if their are rather 
delicate differences in interactions. However, in this 
context, if a physical phenomenon can be characterized 
within a universal concept, such that the chemical de¬ 
tails only contribute to a pre-factor, then the correct 
physics can be captured with a rather simple generic 
model. The use of generic schemes has several advan¬ 
tages; 1) the parameter space is not restricted to a spe¬ 
cific system, unlike the all-atom simulations, and a broad 
range of systems can be represented within a unified sim¬ 
ulation protocol, 2) the time scale of simulations are not 
of concern, and 3) because of the absence of any com¬ 


peting energy scales, one does not need an advanced 
molecular dynamics scheme. In this context, we have re¬ 
cently shown that the complexity of smart polymers can 
be captured within a generic model [lj. Using a simple 
model we could quantitatively capture the reentrant coil- 
globule-coil scenario of PNIPAm and PAPOMe [HI 
in aqueous methanol mixtures. Our analysis suggested 
that when two competing and individually good solvents 
are mixed together, because of the preferential binding 
of the better of the two (co)solvents with the polymer, 
it collapses within the intermediate solution composi¬ 
tions. At a low cosolvent concentration, the cosolvent 
molecules can bind to two distinctly far monomers form¬ 
ing bridges and leading to polymer collapse. When the 
concentration of the better cosolvent is increased, they 
decorate the whole polymer and the polymer opens up. 
These results are in good agreement with the simula¬ 
tions incorporating all atom details [l4j and experiments 
Q. This coil-globule-coil scenario is a generic effect and 
many polymers are expected to exhibit similar behavior 
as long as one of the solvents is significantly better than 
the other. Thus the behavior is not strictly restricted to 
the so called smart polymers exhibiting a LOST. In Table 
[Qwe present a list of polymers that show co-non-solvency 
effect. It is interesting to note that well known stan¬ 
dard polymers, such as poly(ethylene oxide) (PEO) and 
polystyrene, also show co-no-solvency [la 17]. Another 
example may include polymeric semiconductors [22 : ] that 
show anomalous viscosity with solvent composition, sug¬ 
gesting a change in polymer conformation. 

One of the most intriguing aspects of the co-non¬ 
solvency effect is that the solvent quality becomes in¬ 
creasingly better by the addition of the better cosolvent. 
Thus the polymer collapses in a good solvent, making the 
solvent quality decoupled from the polymer conforma¬ 
tion. This is very striking and against the conventional 
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TABLE I: A table listing various polymer systems that show co-non-solvency effect when solvated in their respective mixture 
of solvents. 


Polymer (p) 

Solvent (s) 

Cosolvent (c) 

Poly(N-isopropylacrylamide) (PNIPAm) [6-_9, ^8] 

Water 

Methanol, Tetrahydrofuran, 
or 1,4-dioxane 

Poly(acryloyl-L-proline methyl ester) (PAPOMe) [1A] 

Water 

Methanol 

Poly(ethylene oxide) (PEO) [16] 

Water 

N, N-dimet hy lformamide 

Polystyrene [17] 

N,N-dimethylformamide 

Cyclohexane 

Poly(vinyl alcohol) [HI] 

Poly(2-(methacryloyloxy)ethylphosphorylcholine) 

Water 

Dimethyl sulfoxide 

(PMPC) [20, .21] 

Water 

Methanol, Ethanol, or Iso-propanol 


view on polymer solutions. While the atomistic simula¬ 
tions [1J] clearly demonstrate that the solvent becomes 
increasingly better by the addition of better cosolvent, 
it is still difficult to identify the preferred local coordi¬ 
nation and especially the bridging. In contrast, generic 
simulations give a clear microscopic understanding of this 
complex phenomenon within a simple simulation proto¬ 
col Q. 

Complementary to that computer simulations give a 
good microscopic picture of the polymer collapse tran¬ 
sitions, it is also advantageous to devise a general ana¬ 
lytical theory consistent with the findings known from 
computer simulations and/or experiments. Moreover, 
because of the complexity of the system interactions, 
this discrete particle based phenomenon can not be ex¬ 
plained using a Flory-Huggins type mean-field picture. 
Instead, these complex conformational transitions can be 
explained within a Langmuir-like thermodynamic treat¬ 
ment of competitive displacement of different solvent 
components onto the polymer DJ. 

In this work, we revisit the co-non-solvency effect of 
smart polymers in the mixtures of solvents. We extend 
the analysis of our previous work T| to better under¬ 
stand the microscopic picture of co-non-solvency. We will 
present an in-depth argument to show that the mean-field 
theory is highly unsuitable for these systems and the con¬ 
ceptual need of a discrete particle-based theory. We also 
propose a phase diagram to identify the conformational 
states of smart polymers in various bulk solutions and 
with the change of cosolvent concentrations. 

The remainder of the paper is organized as follows: 
in Section eh the generic molecular dynamics simulation 
details are presented. Results and the theoretical argu¬ 
ments are presented in Section fllll and we finally present 
our conclusions in Section ED 


II. MODEL AND METHODS 

We start by briefly describing the details of the generic 
molecular simulations. A similar model has been used in 
our earlier study. A detailed description of model and 


method is presented in Ref. [lj. Here a polymer p is mod¬ 
eled using the well known bead-spring polymer model 
[23l |. In this model, individual monomers of a polymer in¬ 
teract with each other via a repulsive 6-12 Lennard-Jones 
(LJ) potential (WCA potential). Additionally, adjacent 
monomers in a polymer are connected via a finitely ex¬ 
tensible nonlinear elastic potential (FENE). Here p — p 
interaction energy is chosen as e p = l.Oe and the size of 
the monomer is a p = l.Oer. All units are expressed in 
terms of the LJ energy e, the LJ radius a, and the mass 
m of individual particles. This leads to a time unit of 
r = <7y/ m/e. The parameters of the potential are such 
that a reasonably large time step can be chosen, while 
bond crossing remains essentially forbidden. 

A bead-spring polymer is solvated in mixed solutions 
composed of two components also modeled as LJ beads, 
solvent s and cosolvent c, respectively. Since the solvent 
molecules typically are much smaller than the monomers 
of PNIPAm and/or PAPOMe in aqueous methanol, we 
choose the sizes of (co)solvents to be a s / c = 0.5ct. Note 
that because of the reduced size of (co)solvents, the cor¬ 
responding number density within the simulation domain 
should also be adjusted such that the overall pressure re¬ 
mains ~ 40e/er 3 . p — s and p — c interactions are chosen 
such that c is always a better solvent than s. In our ear¬ 
lier study [lj, the default system consisted of a repulsive 
p — s interaction, while p — c interaction was attractive. 
In this work, we generalize this and investigate the effect 
of assymetry in interaction energies, when both p— s and 
p — c are attractive with interactions e ps and e pc , respec¬ 
tively. Additionally, we impose conditions; (1) e ps < e pc , 
(2) 0.5 < £ ps < 1.0 and 0.5 < e pc < 2.5. Temperature is 
set to T = 0.5 e/k;b, where kb is the Boltzmann constant. 
This leads to a relative energy scale e pc — e ps < 3 kbT. 
These values are typically comparable to the interaction 
energy scale for PNIPAm in aqueous methanol. Solvent 
particles always repel each other with a repulsive LJ po¬ 
tential, with Cij = l.Oe. This is a good approximation 
given that the p — c and p — s interactions are dominant, 
which will be discussed at a later stage. 

The cosolvent mole fraction x c is varied from 0 (pure s 
component) to 1 (pure c component). We consider three 
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FIG. 1: (color online) Normalized radius of gyration R g = 
Rg/R g (x c = 0) as a function of cosolvent molar concentration 
x c for three different chain lengths TV;. Results are shown 
for the default system taken from Ref. [1]. The error bars 
are the standard deviations calculated from six independent 
simulations. The lines are drawn to guide the eye. 


different polymer chain lengths TV; = 10, 30 and 100, 
solvated in 2.5 x 10 4 solvent molecules for TV; = 10 & 
TV; = 30 and 10 x 10 4 solvent molecules for TV; = 100, re¬ 
spectively. The equations of motion are integrated using 
a velocity Verlet algorithm with a time step St = 0.005r 
and a damping coefficient T = 1.0 t~ 4 for the Langevin 
thermostat. The initial configurations are equilibrated 
for typically several 10 5 r, depending on the chain length, 
which is at least an order of magnitude larger than the 
relaxation time in the system. After this initial equili¬ 
bration averages are taken over another 10 4 r to obtain 
observables, especially gyration radii T? g , chemical po¬ 
tentials fx p of the polymer and the bridging fractions of 
cosolvents ^>b- 


III. RESULTS AND DISCUSSIONS 
A. Co-non-solvency: A brief overview 


In our previous paper , lj, we provided a possible expla¬ 
nation for the experimentally observed co-non-solvency 
effect of smart polymers in aqueous mixtures In 

Fig. 0 we show the normalized radius of gyration R g = 
R g /R g (x c = 0) as a function of cosolvent molar fraction 
x c . The data is shown for three different TV;. It can be 
appreciated that just by adding a small fraction of the 
better of the two solvents, the polymer collapses into a 
compact globule structure. As discussed in the intro¬ 
duction, this reentrant collapse and swelling transition is 
facilitated by the preferential binding of cosolvent com¬ 
ponents with the polymer. The initial collapse is due 
to the formation of bridges that the cosolvent molecules 
form by binding two monomers that can be distinctly far 
along the backbone of a polymer, while the reopening at 



FIG. 2: (color online) Static structure factor S(q) for a chain 
length TV; = 100 and for two different mole fractions x c . A 
power law of q -1 '" with v = 0.6 shows an extended (good 
solvent) conformation and q -4 supports a compact globule 
structure. For comparison, we have also plotted the analytical 
scattering function of a sphere. 


higher concentrations is due to the increased decoration 
of the polymer by cosolvent molecules. Therefore, we can 
identify two kinds of cosolvents: fraction <f >b of bridging 
cosolvents that bind to two monomers and a fraction <f> 
of cosolvents that are only bound to one monomer. Note 
that other than cf>B and q !>, there are a large fraction of 
free cosolvents that are present in the bulk solution and 
usually are required to maintain solvent equilibrium. 

It is also interesting to observe an inverse system size 
effect in the reopening transition, as observed in Fig. O 
While the initial collapse (x c < 0.1) is reminiscent of 
a first-order-like collapse, the reopening (for x c > 0.5) 
is rather smooth even for longer chain lengths. This is 
contrary to the knowledge of critical phenomena. Thus 
indicating that this transition is not a phase transition 
in a true thermodynamic sense. This aspect will be dis¬ 
cussed at a later stage of this manuscript. Furthermore, 
the cosolvent driven first-order-like collapse for x c < 0.1 
is reminiscent of the temperature induced first order tran¬ 
sition in PEO [24| . 

To further quantify the collapse facilitated by bridging 
cosolvents, we calculate the static structure factor S(q) 

mm, 



(i) 


In Fig. Owe present S(q) for TV; = 100. As expected, 
a power law well approximated by q -1 / 0 - 6 is observed 
for x c = 0 (pure solvent), a signature characteristics of 
an extended coil structure. For x c = 0.1, the polymer 
collapses into a compact globule, with ~ 60% decrease in 
R g with respect to its original extended R g at x c = 0 (see 
Fig-0), as shown by a prominent scaling law q 4 in Fig. [2] 
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[25j . Note that R g for x c = 0.1 is slightly larger than the 
equivalent R g when a polymer collapses because of pure 
depiction effects. This is due to the fact that a collapsed 
polymer also contains interstitial bridging cosolvent and 
their sizes contribute towards a slightly larger compact 
globule. 

Preferentiability is required for the observation of the 
co-non-solvency effect. Therefore, conformations of poly¬ 
mers in mixed solvents are intimately linked to the asym¬ 
metry in p — c andp— s interactions £ pc —£ ps - In Fig.[3]we 
present a unified picture of polymer conformation with 
changing £ pc — £ ps at different x c . It can be appreciated 
that, for £ ps = £ pc , co-non-solvency is not observed and 
the polymer remains in a coil conformation. Only when 
£ P c ~ £ps > 0.25Kb? 1 , does the polymer exhibit a coil- 
globule-coil-like scenario. More interestingly, the larger 
the difference £ pc — £ ps the smoother the re-opening tran¬ 
sition at larger x c values. This is not surprising given that 
for a stronger £ pc , </>b has stronger binding, thus leading 
to a more stable semi-collapsed conformation. Thermo¬ 
dynamically, the energy density in the solvation shell can 
be increased by increasing the p — c interaction strength. 
Increasing energy density by a factor of two will approx¬ 
imately act in a similar manner as that of a polymer of 
twice Np Interestingly enough, increased £ pc and/or N( 
has the same effect on the overall polymer conformation. 

We also want to point out that, for £ pc » £ ps , poly¬ 
mer collapse will occur close to x c —> 0. A more promi¬ 
nent representation of the region 0.0 < x c < 0.1 will 
require fine grids and systematic scanning of the concen¬ 
trations within the range 0 < x c < 0.1. Here, however, 
while the initial collapse is always first order like, the 
re-opening has much stronger dependence on £ pc — £ ps . 
Therefore, we abstain from presenting any more details 
for 0 < x c < 0.1. 

The most interesting aspect of this reentrant transi¬ 
tion is that even when the solvent quality becomes bet¬ 
ter and better by the addition of the better good solvent, 
the polymer collapses in good solvent. This makes the 
polymer conformation decoupled from the solvent qual¬ 
ity and only dictated by the preferential coordination of 
cosolvent with polymer. This particle based phenomenon 
can not be explained within a mean-field type approach. 

Before describing an analytical theory, we briefly want 
to comment on the suitability of the generic simulation 
protocol to study the complexity of smart polymers. One 
important aspect of smart polymers, such as PNIPAm, 
PVCL, and PAPOMe, is their thermal responsiveness. 
These polymers remain in a coil configuration at low tem¬ 
peratures, while collapsing into a compact globule at high 
temperatures, thus presenting a lower critical solution 
temperature (LOST). In this context, it is worth men¬ 
tioning that the generic schemes do not present LOST 
and co-non-solvency can be studied at one fixed tem¬ 
perature over full range of x c . Moreover, the co-non¬ 
solvency effect is not necessarily restricted to the “so 
called” smart polymers exhibiting LCST. Therefore, a 
broad range of polymers are expected to show a simi¬ 


lar reentrant scenario, as long as they are dissolved in a 
mixture of competing good solvents. A list of the pos¬ 
sible polymer systems that show co-non-solvency is pre¬ 
sented in Table |U Standard polymers, such as PEO in 
aqueous DMF [Hi and polystyrene in a mixture of N,N- 
dimethylformamide (DMF) and cyclohexane [17], also 
show co-non-solvency [n|. Another example includes 
polymeric semiconductors. It has been observed that the 
solution viscosity of polymeric semiconductors can dis¬ 
play a non-monotonic dependency on cosolvent concen¬ 
tration, indicating conformational change [22]. There¬ 
fore, we speculate that many more polymers, such as 
polycarbonate or polypropylene, may also exhibit a sim¬ 
ilar reentrant transition in appropriate competing good 
solvents. 

Another aspect is that the NIPAm monomer has a 
hydrophilic part and two hydrophobic parts, as shown 
in the schematic Fig. [I] It is generally believed that 
methanol molecules bind to the hydrophilic part and 
thus push away water molecules towards the hydropho¬ 
bic part, leading to polymer collapse. In a generic simu¬ 
lation, however, the monomer is represented by a sphere, 
thus eliminating any effects due to hydro(phob/phil)icity 
within the model. Our simulations (Ij, [3] suggests that 
the only dominant interaction is the preferential coor¬ 
dination of methanol around the NIPAm monomer (see 
Fig. 3 of Ref. 0 )- The collapse is initiated by bridging 
and not by any hydrophobic effects that may occur due to 
solvent interaction with the alkane backbone. Chemical 
details do not play any role in describing this reentrant 
transition. In this context, tuning specific (co)solvent- 
polymer interactions, a whole new class of heteropoly¬ 
mers can also exhibit co-non-solvency and related phe¬ 
nomena (26| . 


B. Co-non-solvency: A simple analytical approach 

1. Why mean-field theory is inappropriate to describe 
co-non-solvency 

When a polymer with chain length A) at volume frac¬ 
tion cj) p is dissolved in a mixture of two components s and 
c, respectively, the standard Florv-Huggins energy Jth 
of polymer solutions reads [23,123] , 

— -f; = tt In <t>p + x c (1 - </> p ) In [x c {l- <j) p )\ 
kb? iv i 

+ (1 - x c ) (1 - (j) p ) In [(1 - x c ) (1 - (f> p )] 

+ X PS <M 1 - x c ) (1 - </> P ) 

A X P c l l ) P x c (1 /’p) 

+ XscXc (1 - X c ) (1 - (j>p) 2 . (2) 

Here, the first three terms represent the entropy of mixing 
and the last three terms deal with interactions between 
different components i and j via Xij- Expanding Eq.[2]to 
the second order gives a direct measure of the excluded 
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FIG. 3: (color online) A sketch of the phase diagram showing the change in normalized radius of gyration R g = R g /R g (x c = 0) 
with the varying cosolvent molar concentration x c and the relative interaction strengths e pc — £ ps . Results are shown for chain 
length Ni = 30. In the contour plots, the area bound by the black curve represents maximum collapse with R g < 0.6. The red 
contour curve represents the boundary when polymer goes from globule-coil (or vice-versa) by either changing x c at a constant 
E pc — e ps or by changing E pc — £ ps at constant x c - The region outside the dashed green curve shows maximum extension of the 
chain with R g > 1.0. Prominent kink in the black contour curve is due to the error bar associated with that data points. 



(b) 

phobic 


philic 
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FIG. 4: (color online) Part (a) shows simulation snapshot 
representing a monomer of PNIPAm. Hydrogen atoms is ren¬ 
dered in steel, green spheres are Carbon atoms, blue sphere 
is Nitrogen, and the Oxygen is rendered in red. Part (b) 
represents a schematic representation of NIPAm monomer. 


volume V of the polymer [13, HU ; 

V — 1 2(1 x c ) Xps ‘^'EcX.pc T 2x c (1 x c ) Xsc(3) 

where Xps and Xpc are the Flory-Huggins interaction pa¬ 
rameters between p — s and p —c, respectively. The factor 
Xsc is the parameter of s — c interaction. When both sol¬ 
vent and cosolvent are good solvents, Xps < 1/2 and 
X pc <1/2 @]. Using the first two terms of Eq. [3] 


we find a linear variation of V with x c for the cases 
of non-interacting s and c (i.e. Xsc = 0), as shown 
by the blue line in Fig. [5] It is also clear from Fig. [5] 
that only when Xsc < 0 can V become negative, opening 
the possibility for the coil-to-globule-to-coil conformation 
changes typical of co-non-solvency. It has been noticed 
early 0] that for common solvent mixtures where co-non¬ 
solvency effects are observed, such as water-alcohol mix¬ 
tures, Xsc > 0, thus precluding any explanation based on 
a mean-field, Flory-Huggins type of analysis. 

Furthermore, within the mean-field picture described 
in Eq. [2l one can get the expression for the shift in chem¬ 
ical potential of polymer p, p for <j> p —» 0, 


Mp (<^p t 0) 


dT FH 


d<j) p 


n —>-0 


= const — x c In x c 


(1 - x c ) In (1 


T (1 %c)Xps “F %cXpc 

2x c (1 x c ) Xsc- 


X c ) 


(4) 


In Fig. [Gl we present a schematic representation of p, p as 
expected from Flory-Huggins picture described in Eq. [I] 
Here Jt p (x c = 1) < p. p {x c = 0) because alcohol is a better 
solvent compared to water. Consistently with the behav¬ 
ior of V presented in Fig. [5j p, p for Xsc < 0 displays a 
hump for intermediate mixing ratios where the solvent 
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FIG. 5: (color online) A schematic representation of the poly¬ 
mer excluded volume v as a function of cosolvent mole fraction 
x c . The curve shows that the interaction parameter \' sc be¬ 
tween solvent-cosolvent is a key factor to exhibit a swelling- 
collapse-swelling scenario. When \sc = 0 the polymer re¬ 
mains swollen. 
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FIG. 7: (color online) Kirkwood-Buff integral Gij between 
different solution components as a function of cosolvent molar 
fraction x c . Lines are the polynomial fits to the data that 
are drawn to guide the eye. The data was obtained from 
the semi-grand canonical simulations incorporating all-atom 
details [141 ]. For pure solvent at x c = 0.0 and pure cosolvent at 
x c = 1.0, individual coordinations G pc and G ps are undefined, 
respectively. 

p — c interactions to observe this reentrant transition. 
Moreover, if the mean-field theory is sufficient to under¬ 
stand this reentrant coil-globule-coil transition then the 
analysis of the bulk solution property, that can easily be 
calculated using molecular simulations, should also show 
a preferred s — c coordination over p — s or p — c coordi¬ 
nation. 

A quantity that best describes the relative intermolec- 
ular affinity and/or the interaction strength is the fluc¬ 
tuation theory of Kirkwood and Buff (KB) [29j. KB the¬ 
ory connects the pair distribution function to thermody¬ 
namic properties of solutions using the “so called” KB 
integrals; 


FIG. 6: (color online) A schematic representation of the shift 
in chemical potential JT p as a function of cosolvent mole frac¬ 
tion Xc- 


quality goes from good to poor to good again (see red 
curve in Fig. El) . However, in our simulations [1], |l4[ not 
only is Xsc = 0, but we also measure a chemical poten¬ 
tial trend similar to the black schematic curve in Fig. [6] 
Thus the solvent quality remains good in the whole com¬ 
position range and, in-fact, it even becomes increasingly 
better by the addition of the cosolvent. By the analysis 
of Eq. U it can be seen that a similar trend as the black 
curve of Fig. [S]can be obtained from the mean-field pic¬ 
ture when Xsc >> 0. However, this can only be obtained 
at the nonrealistic cost of driving the system towards 
solvent phase separation. This further confirms the in¬ 
capability of mean-field theory to capture the reentrant 
co-non-solvency effect in polymeric systems. 

The mean-field picture also suggests that the strength 
of s — c interaction should be dominant over p — s and 


Gij = 4tt / [g ij(r) - 1] r 2 dr, (5) 

Jo 

where g ij(r) is the pair distribution function. In Fig. [7] 
we summarize Gij between different solvent components. 
It can be appreciated that, for 0.1 < x c < 0.5, p — c co¬ 
ordination is at-least an order of magnitude larger than 
the values between the solvent components in the 
bulk solution, suggesting that the fraction of cosolvent 
molecules in close contact with the chain is always much 
larger than its natural, mean-field proportions in the bulk 
solution. This is contrary to what is know from the analy¬ 
sis based on the mean-field theory presented above. Fur¬ 
thermore, the shift in chemical potential fi can be es¬ 
timated from the KB theory. If a polymer p at dilute 
concentration is solvated in a mixture of solvent s and 
cosolvent c, p p can be calculated using [30], 

^Ee\ _ G ps — G pc /g\ 

d Pc )p T 1 — p c (G cs — G cc ) ’ K) 

where ]I p = p p / k-bT , and p c is the cosolvent number den- 
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FIG. 8: (color online) Shift in chemical potential of a single 
monomer JL as a function of methanol mole fraction x c . Jl p 
is calculated by integrating the Eq. [6j Gij are taken from 
Fig. [7] Dashed line is plotted according to the Eq. [13] 

sity. The change in ~p p is shown in Fig. [5] Data clearly 
show the trend known for the case when \ sc = 0, suggest¬ 
ing that the solvent quality becomes better and better by 
addition of the better (co)solvent. This decoupling be¬ 
tween solvent quality and the polymer conformation is 
contrary to the conventional understanding from mean- 
field predictions. Additionally, this conformational tran¬ 
sition of polymer in mixtures of competing good solvents 
is not a phase transition in true thermodynamics sense 
and is only dictated by the preferential adsorption of one 
of the (co)solvents. Therefore, an analytical description 
is needed that can incorporate the concept of competitive 
adsorption by taking into account the strong deviations 
of local concentration from mean-field values. We will 
address this in the following section. 

2. Competitive adsorption of the cosolvent as a model for 
co-non-solvency 

We have recently proposed that the polymer collapse in 
co-non-solvency phenomena can simply be understood as 
the result of the attractive interactions induced by cosol¬ 
vent molecules that form a bridge between two monomers 
[l]. Conformational collapse, at low x c , is thus in¬ 
duced by the increase of such bridges, while polymer 
swelling at larger cosolvent fractions is due to the pro¬ 
gressive replacement of these bridges by single site cosol¬ 
vent molecules that are attached to one monomer only. 
Thus, for high enough x c , these non-bridging cosolvent 
molecules eventually decorate the whole chain backbone 
to facilitate the reopening. In Fig. [9] we show </>b, the 
fraction of backbone sites participating in bridge forma¬ 
tion, as a function of x c . Note that cf and <f >b are cosolvent 
molecules that are directly in contact with the monomers 
at a distance 2 1 / 6 o- pc ~ 0.84er. It can be appreciated 


FIG. 9: (color online) Bridging fraction of cosolvents </>b as 
a function of cosolvent mole fraction x c for Ni = 100. The 
data corresponding to red o is the direct calculation of q !>b 
from the simulation trajectory. The prediction of analytical 
theory from Eq. [8] is plotted for two different parameters. 
Here £ = 0.05 corresponds to the translational entropic term 
corrected with a loop contribution, as observed earlier [jj. 

that, within the range 0.1 < x c < 0.4, <f >b obtained 
from the numerical simulations shows a distinct hump 
that is consistent with the range of x c when the poly¬ 
mer collapses into a (compact) globule and then gradu¬ 
ally begins to expand. To devise a theoretical formula¬ 
tion, we view the polymer as a substrate with N sites 
exposed to the bulk solution, of which n s sites are occu¬ 
pied by s (solvent) molecules, n c sites by non-bridging c 
(co-solvent) molecules and 2 n B sites by bridging c (co¬ 
solvent) molecules, with N = n s +n c +2n c B . The observed 
sequence of collapse and re-swelling of the polymer cor¬ 
responds to a fast growth of n c B as x c increases, followed 
by a displacement of n B by n c for larger x c values. Such 
a sequence is typical for competitive displacement in ad¬ 
sorption phenomena [31] . Our results from numerical 
simulations for n c B and n c , or alternatively for the frac¬ 
tions 4>b = n c B /N and <f> = n c /N, are very well described 
by a competitive adsorption model with the following as¬ 
sociated free energy density of adsorption for non-bridges 
and bridges, 

—In (cf) + (cj) B In (2(f B ) 

k b T 

+ (1 — (f - 2<£b) In (1 -<j> — 2(j) B ) 

- £4> — £ B ^) B -^=(<£ + </>b)) (7) 

k b 2 

with pi = K B Tln(x c ) being the chemical potential of the 
cosolvent in the bulk solvent mixture and the adsorption 
energies £ and £ B measure the excess affinities of individ¬ 
ual non-bridging and bridging cosolvent molecules to the 
chain monomers. The first three terms in Eq. [7] express 
entropic contributions of the adsorbed bridges and non¬ 
bridges to the energy densities, while the two following 






terms measures contact energies between the cosolvents 
bridges and non-bridges with the polymer backbone. The 
unusual pre-factor £ is, as discussed later, a consequence 
of assuming a logarithmic form for the dependence of the 
energy required to make a bridge on the average density 
of existing bridges. This is the case for instance [lj, if one 
assumes that in order to make a new bridge at density 
0 b, the chain needs to make a loop of length £ = 1/0b, 
with associated penalty ~ \og£ ~ log(l/0B). 

Minimization of Eq. [T] with respect to 0 b and 0 leads 
to the implicit equation for the bridge density 4>b(x c ), 

16^ x c = (i — 20 b) 

± . ( 8 ) 

with x* = exp (—£) and x** = exp (—£b + 21 n 2 e — £) 
are the characteristic concentrations related to the ad¬ 
sorption energies £ and £b for non-bridges and bridges. 
Fig. [5] shows that this expression describes very well our 
experimental results, with £ = 0.05. 

Eq. [ 8 ] can equivalently be derived by considering the 
two pseudo chemical reactions, 

cosolvent + empty site non — bridge 
cosolvent + 2 empty site ^ £ bridge. (9) 

A schematic representation of this reaction is presented 
in Fig. [10] When the solvent and cosolvent interactions 
with the polymer backbone empty sites are described as 
pseudo reactions, a cosolvent molecule reacts with one 
empty adsorption site to form one adsorbed non-bridge, 
while it reacts with two empty sites to make £ bridges. 
The associated equilibrium standard mass-action laws 
can thus be written as 

Xc _ <t> 

x* 1 - 0 - 20 B 

-?£- = _ ^ _ (10) 

4®S* (1-0-2</> b ) 2 ’ 1 ; 

with equilibrium reaction constants l/x* and 1/x**. 
Note that the reaction equilibrium concentration x** has 
been, for mathematical convenience, defined up to a fac¬ 
tor four. Solving the mass-action laws for 0b gives Eq. [8] 
In this pseudo-chemical language, the factor £ describing 
the effective number of bridges formed by the interaction 
between one cosolvent molecule and the two empty sites 
of the backbone appears as a consequence of assuming a 
power-law dependence for the equilibrium constant of the 
pseudo-chemical reaction. Note that the actual shape of 
Eq. [5] is quite sensitive to the value of £. In particular, 
the choice £ = 1 , corresponding to a standard chemical 
reaction between free species in solution, leads to a pre¬ 
diction that can not describe our data (see red curve in 
Fig. EJ. 
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(b) } + « = J, 


+ (5 = 

••• 

FIG. 10: (color online) A schematic representation of the 
chemical reaction described in Eq. El Part (a) describes a 
typical polymer conformation decorated by non-bridging and 
bridging cosolvent molecules. While part (b) shows that a 
polymer segment and a cosolvent forms a single adsorbed non¬ 
bridging cosolvent, part (b) represents two segments and a 
cosolvent makes a bridge (or a bridging cosolvent). 

In a previous work (lj, we argued that a value of 
£ = 0.05 can be understood by considering loop contri¬ 
butions to the cost of making a bridge. When a pure con¬ 
figurational cost for distributing the bridges amongst the 
possible occupation sites is combined with the entropic 
cost of loop formation, one can write £ = 2 — m. Here 
the critical exponent m can be estimated within a simple 
scaling argument. In this context, one can characterize 
the loop formation by a partition function of vanishing 
end-to-end distance R e — > 0 [28j ] . 

Z Nl (R e ^0)ocq N ‘N l a - 2 , ( 11 ) 

and the partition function at finite R e is given by, 

Z^R^cxq^N^- 1 . (12) 

Here 1/q is the critical fugacity and the universal expo¬ 
nent a = 0.2 HI. From these two cases one can esti¬ 
mate the free energy barrier to form a loop of length t as 
AT(£) = hikbT ln(^), with m = 7 — a + 1 being the crit¬ 
ical exponent Hi. Although this gives m — 1.95 for loop 
formation in self-avoiding walks, in excellent agreement 
with our findings, it is worth pointing to the fact that 
our simple analytical description does not address other 
possible contributions to bridge formation, such as coop¬ 
erative or other non-trivial entropic effects that might be 
determinant in the dense chain globule. 

This selective adsorption model provides also for an 
analytical prediction of the shift in the chemical potential 
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fi p as a function of x c , 


Pp 

kb T 


const + (2 — £) 4 >b 


In 


1 + </>B 


1—C/2 




Fig. [8] shows a comparison between predictions from 
Eq. [13] and the values of the chemical potential obtained 
from Eq. [6] A very good agreement is obtained by simply 
inserting into Eq. [13] the values for £ and concentrations 
obtained from the fit of the bridging fraction, further 
confirming the consistency and validity of our approach. 


IV. CONCLUSION 

We present a comprehensive analysis of the co-non¬ 
solvency effect of smart polymers in a mixture of good 
solvents. Our results suggest that co-non-solvency is a 
generic effect that is not restricted to any specific chemi¬ 
cal systems. Though co-non-solvency has been associated 
with polymers like PNIPAm, PVL, and/or PAPOMe 
0, ii, Ea El, polymers such as PEO and polystyrene 
also exhibit co-non-solvency HEHil- Furthermore, this 
reentrant transition is dictated by the preferential co¬ 
ordination of one of the cosolvents with the polymers. 


More interestingly, even when the chain collapses, the 
solvent quality becomes increasingly better. This makes 
the solvent quality disconnected from the conformation 
of the macromolecules. This discrete particle-based phe¬ 
nomenon can not be explained within a mean-field the¬ 
ory. Instead it can be explained using a thermody¬ 
namic treatment of a simple selective adsorption picture 
and a generic simulation protocol. Therefore, this work 
presents a unified theoretical and computational frame¬ 
work, which can pave the way for a more generic under¬ 
standing of polymeric solubility in mixtures of solvents. 
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